Evidence for orbital ordering in LaCoO^ 
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We present powder and single crystal X-ray diffraction data as evidence for a monoclinic distortion 
in the low spin {S = 0) and intermediate spin state (5* = 1) of LaCoO^. The alternation of short 
and long bonds in the ab plane indicates the presence of Cg orbital ordering induced by a cooperative 
Jahn- Teller distortion. We observe an increase of the Jahn- Teller distortion with temperature in 
agreement with a thermally activated behavior of the Co^^ ions from a low-spin ground state to an 
intermediate-spin excited state. 
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The study of orbital degrees of freedom in transition 
metal (TM) oxides has gained prominent interest. Novel 
techniques such as Resonant X-ray Scattering and X-Ray 
Absorption Spectroscopy give direct information about 
orbital occupancy. It is realized that the orbital moments 
are as important as the spin moments to understand the 
electronic properties. Prominent examples in which the 
orbital degrees of freedom determine the electronicprop- 
erties are the metal-insulator transitions in V2O3 ll| and 
doped LaMnO'i 2\. Here, spin- and orbital- induced tran- 
sitions are intimately related. In the perovkites, metal- 
licity and orbital ordering are mutually exclusive because 
of the large Jahn- Teller (JT) splitting of the Cg orbitals 
0- Furthermore, novel ground states can be obtained 
such as an orbital liquid 0, 0] . 

Considerable interest in the perovskite-type LaCoO^ 
originates from the puzzling nature of two transitions 
in this compound and the vicinity to a metal-insulator 
transition. The ground state of LaCoO^ is a nonmag- 
netic insulator and there is no long range magnetic order 
at all temperatures. At low temperatures, the magnetic 
susceptibility increases exponentially with temperature 
exhibiting a maximum near lOOK. At higher tempera- 
tures, a second anomaly is observed around 500K which 
is accompanied by a semiconductor to metal transition. 
The maximum at lOOK was ascribed initially to a change 
of the spin state in the Co^+ ions i.e. a transition from 



a low-spin (LS) nonmagnetic ground state 
to a high-spin (HS) state (tLel^S = 2^ " 
the more recent literature Il3, nil lia, [lA |l 



t^^S = 0) 
1- In 

new sce- 
narios involving an intermediate-spin state (IS) (^Ig^gj 
5 = 1) have been proposed. Using LDA-I-U calculations, 
Korotin et al. [la| proposed the stabilization of the IS 
state due to the large hybridization between the Co-Cg 
and 0-2p levels. Due to the partially filled Cg level, the 
IS state is JT active. The degeneracy of the eg orbitals 
of Co'^"'" ions in the LS state is expected to be lifted in 
the IS state by a JT distortion. 

All structural studies, based on powder X-ray and 
neutron diffraction experiments are consistently inter- 
preted in rhombohedral R3c symmetry and no struc- 
tural transitions are reported in the temperature interval 



4.2K - 1248K E [H H [13 . A cooperative JT distor- 
tion is incompatible with this space group. The rhom- 
bohedral distortion of the parent cubic perovskite struc- 
ture consists of a deformation along the body diagonal, 
and preserves only one Co — O distance. This triggered 
our present study to detect a coherent JT distortion and 
to investigate why recent high quality structural studies 
failed to observe such state. As the JT distortion caused 
by eg orbital ordering are usually quite large, there must 
be reasons why the distortion went unobserved. Further- 
more, the LS-IS transition is robust and observed both in 
powder and single crystal samples. A lattice distortion 
due to the JT effect in the IS state has been suggested 
by infrared spectroscopy measurements 0| and recently 
by magnetic susceptibility and thermal expansion mea- 
surements 20]. 

It is in certain cases exceedingly difficult to identify 
the correct crystal symmetry. In perovskites such as 
LaAInOs or YVO3 the JT distortion involves alternating 
long and short TM-0 distances without always lowering 
the symmetry of the lattice [ij. Consequently, refine- 
ment of powder diffraction data are not very sensitive to 
several symmetry modifications and single crystal diffrac- 
tion is prerequisite to determine the crystal symmetry. 
LaCoOs has been consistently analyzed with the rhom- 
bohedral space group R3c. This rhombohedral symme- 
try involves an alternating rotation of the corner sharing 
CoOq octahedra along all three crystallographic axes of 
the undistorted, cubic perovskite parent structure. The 
corner sharing oxygen octahedra can be described by one 
oxygen position. This symmetry allows only one Co — O 
distance and thus rules out a coherent JT distortion. 
Such distortion, while preserving the alternating rotation 
of the CoOq octahedra, lowers the symmetry to mono- 
clinic subgroups of R3c. The monoclinic space group 
12 /a is a subgroup of i?3c with the highest symme- 
try that accommodates the rotations of the octahedra 
as well as a coherent JT distortion. The symmetry re- 
duction to 12 /a space group is caused by rotations of 
unequal magnitude and/or a differentiation of Co — O 
bond lengths, induced by a coherent JT distortion. In 
this symmetry two oxygen positions describe the oxygen 
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framework. A shift of the second oxygen position with 
2^02 = yo2 corresponds with a differentiation of Co — O 
bond lengths, whereas a movement perpendicular to this 
constraint results in a coherent rotation of the octahedra. 

We provide strong evidence for such monoclinic distor- 
tion in the LS /IS state using high resolution synchrotron 
single crystal and powder diffraction experiments. The 
monoclinic distortion is ascribed to a JT distortion that 
signals the long range ordering of the Cg orbitals. Our 
results are in agreement with the LS /IS state scenario in 
which the occupation of the IS states of the Co^^ ions is 
thermally activated from the ground LS state ^(\ . 

Single crystals of LaCoO^ were grown by the floating- 
zone technique in an image furnace. High resolution X- 
ray diffraction measurements were performed on a triple- 
axis high energy diffractometer with a photon energy of 
115 keV at the beamhne IIID-C, BESSRC CAT, Ar- 
gonne National Laboratory, us — x scans were collected on 
a single crystal of approximate dimension 1mm. Full data 
sets were measured on a small single crystal fragment of 
approximate dimension 0.1mm using Mo — radiation. 
The measurements were collected using a Bruker SMART 
CCD camera and face indexed absorption corrections 
were performed using associated software packeges^. 
The integration of the frames was carried out using the 
SAINT |2J| software incorporating the absorption cor- 
rection program SADABS. The separation of the twins 
was done using the program GEMINI ji^l- Structural 
determination and refinement was carried out using the 
SHELXTL 13 package. 

We measured the temperature dependence of the mag- 
netization of single crystal LaCoO^. The magnetization 
can be fitted well with a non magnetic ground state 
and a thermally activated S = 1 state (^IgSg) with an 
activation en ergy of A ~ 180K, as shown recently by 
Zobel et al. |2flj. The stoichiometry of the crystal can 
be derived from the very small Curie-term, indicating 
less than 1% S' = 1 impurities. The impurity contribu- 
tion gives rise to the deep minimum at 35K in the mag- 
netic susceptibility. The raising of this minimum in the 
powder samples is attributed to Co^^ localized magnetic 
moments associated with the reconstruction of the sur- 
face 0- The volume of this surface is larger for powder 
samples than for single crystals. Our single crystal mea- 
surement show a minimum susceptibility which is ^ 2 
times smaller than in the powder samples [l^ . 

A full structure determination has been made by sin- 
gle crystal X-ray diffraction. However the crystals are 
highly twinned as is common both for orthorhombic and 
rhombohedral perovskites. Twinning can complicate the 
symmetry determination considerably. 

The rhombohedral perovskites are obtained by a com- 
pression of the parent cubic cell along one of the body 
diagonals. The twinning is a result of the four body di- 
agonals of which only one unique axis is chosen. The 
twin law relates four twin fractions by 90° rotation about 



the cubic main axes. The diffraction pattern of our sam- 
ples typically exhibited two or three rhombohedral twins. 
The twinning "masks" the reflection conditions associ- 
ated with the lowering of the symmetry from i?3c to 12/ a, 
especially the c-glide plane and the three-fold axis. First, 
the violation of the c-glide plane will give rise to extra 
reflections. However, these reflections are also generated 
by the rhombohedral twinning. Therefore, their pres- 
ence cannot be used as evidence for monoclinic distor- 
tion. Second, while the disappearance of the three-fold 
symmetry doesn't generate new reflections, the resulting 
twinning mimics Ric symmetry. A small monoclinic dis- 
tortion can give rise to such pseudo-merohedral twinning, 
which can be easily overlooked in a standard analysis due 
to insufficient angular resolution. Therefore, high angu- 
lar resolution single crystal X-ray diffraction experiments 
are needed to determine the symmetry. 

We verified the presence of twinning by performing lo — 
X scans with high angular resolution (0.01 degrees/step) 
using synchrotron radiation. Such scans measure a cross 
section of reciprocal space on the Ewald sphere perpen- 
dicular to the scattering vector. In this way we measure 
the fine structure of a Bragg peak caused by twinning. 
Fig. 1 shows such a profile of a (10 2 8) reflection [2^ in 
anuj ~x scan measured at 60K. The two big spots which 
split by 0.15(1) degrees in w, correspond to two twins 
generated by lowering the symmetry from cubic to rhom- 
bohedral. Each of these rhombohedral twins shows a fur- 
ther splitting of about 0.04(1) degrees in lo. This smaller 
splitting is best observed at 20K and 60K but also the 
measured rhombohedral peaks at 120K and room tem- 
perature show clear shoulders. We attribute this smaller 
splitting to the twinning caused by the decrease in sym- 
metry from rhombohedral to monoclinic. The same type 
of profiles are observed for all the peaks we have measured 
on several crystals. Such splitting indicates a distortion 
from the rhombohedral unit cell parameters. The ob- 
served splitting indicates that the monoclinic symmetry 
is maintained in the temperature interval 20K - 300K. 

After establishing the symmetry of the structure, the 
atomic positions can be determined from standard single 
crystal diffraction experiments. Full data sets were mea- 
sured at different temperatures (90K, 120K, 140K, 160K, 
200K, 250K and 295K) on a Bruker APEX diffractome- 
ter, using Mo — Ka radiation. The rhombohedral twin- 
ning was clearly observable and the twins were analyzed 
separately. However, the monoclinic twinning could not 
be observed on this diffractometer due to insufficient an- 
gular resolution. The presence of the monoclinic twins 
was observed on a Enraf-Nonius CAD4 diffractometer 
equipped with a point detector which yields a better an- 
gular resolution. Again, these measurements indicate the 
splitting of the rhombohedral peaks as evidence for mon- 
oclinic distortion. However, for structural refinement, 
the monoclinic twin fractions could not be separated. 
Therefore the unit cell parameters derived from the sin- 
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TABLE L Lattice parameters for LaCoOs 



-5.8 -5.7 

m (degrees) 



FIG. 1: The profile of a (10 2 8) reflection |2^ in an - x 
scan at 60K. The two rhombohedral fractions positioned at 
wi ~ —5.82° and TO2 ~ —5.69° split further into three peaks 
representing the monoclinic twins and are labelled as a, b, c 
and d, e, f respectively. 



gle crystal measurements are averaged at their rhombo- 
hedral values (Table 1). The intensities of the reflections 
were obtained from the largest rhombohedral twin frac- 
tion which includes the integrated intensity of the mono- 
clinic twin fractions together. Nevertheless, precise posi- 
tions of the atoms can be obtained by making use in the 
refinements of the monoclinic twin relation. 

The pseudo-merohedral twinning in LaCoO^ results 
from the loss of three-fold symmetry during the crys- 
tal growth or a structural phase transition at elevated 
temperature. Therefore, the twinning corresponds to the 
three possibilities of converting the rhombohedral axes in 
i?3c symmetry into the c-axis in monoclinic symmetry. 
There is no physical reason for one twin to be favored over 
the others and equal volume fractions can be expected. 
All crystals that we investigated, exhibited such a mono- 
clinic twinning. The twin law relates three twin fractions 
by 120° rotation about the cubic [111] direction. If the 
twin fractions are equal, the reflections (h,k,l) and 
cyclic permutations have equal intensities and simulate 
the existence of a three-fold axis. Thus, it is not sur- 
prising that the quality of the refinements in both i?3c 
and 12 /a space groups are essentially the same. How- 
ever, by introducing this twin model, the quality of the 
refinements in 12 /a space group was improved , the re- 
sulting i?l-values falling for a typical data set (140K) 
from 0.0445 to 0.0433. The refined values of the twin 
fractions ranges between 30% and 40%. The atomic co- 
ordinates for a typical data set (200K) and the selected 
bond lengths are displayed in Table II and Table III re- 
spectively. 

Recent high resolution powder neutron diffraction 
measurements were analyzed with i?3c symmetry and 
failed to identify a monoclinic distortion [ijj. We per- 
formed experiments on high quality powder obtained 
from a crushed single crystal using synchrotron radia- 



Volume 



295 5.3695(4) 
250 5.3672(4) 
200 5.3611(3) 
160 5.3578(2) 
140 5.3585(2) 
120 5.3544(2) 
90 5.3469(3) 



5.4331(4) 
5.4329(4) 
5.4316(3) 
5.4291(2) 
5.4316(2) 
5.4294(2) 
5.4227(3) 



7.6395(6) 
7.6367(5) 
7.6318(4) 
7.6288(3) 
7.6304(3) 
7.6257(2) 
7.6168(4) 



90.985(5) 
91.014(4) 
91.056(3) 
91.108(3) 
91.126(2) 
91.148(2) 
91.174(3) 



222.83(4) 
222.65(4) 
222.20(3) 
221.86(2) 
222.04(2) 
221.64(2) 
220.80(2) 



TABLE II; Atomic parameters for LaCoOs at 200K 



Atom Position 


X 


y 


z 


Ui s o 


La 4e 


0.25 


0.25019(13) 





0.00486(11) 


Co 4c 


0.75 


0.25 


0.25 


0.00452(13) 


01 4e 


0.25 


-0.3068 





0.0101(15) 


02 8f 


0.0241(10) 


0.0332(10) 


0.2293(7) 


0.0044(6) 



tion and high angular resolution: 0.001 degrees/step. We 
measured the temperature dependence of the profile of 
the (4 0) reflection 26] as the degeneracy of this peak 
is expected to be lifted by lowering the symmetry from 
rhombohedral i?3c space group to monoclinic 12 j a space 
group. The peak width is resolution limited (0.0025°) up 
to 65K and it starts to broaden gradually as the tempera- 
ture is increased, reaching 0.0036° at 145K. We associate 
the broadening of the Bragg reflection with a gradual 
distortion of the structure. The monoclinic symmetry is 
compatible with a splitting of this reflection. In the left 
side of flg. 2 we show the temperature dependence of the 
splitting of the (4 0) reflection derived from the peak 
proflles. In the inset of flg. 2 we show the profiles of 
the peak at 45K and at 120K. The splitting in derived 
from the lo — x scans on single crystals corresponds with a 
much smaller splitting in — 20 scans in the Bragg geom- 
etry. Thus, due to limitation in angular resolution, the 
monoclinic distortion can go easily undetected in Bragg 
geometry powder X-ray diffraction. 

We notice that powder and single crystal diffraction 
give complementary information. The single crystal anal- 



TABLE III: Selected bond distances in Afor LaCoOs 
Temp(K) Co-02(A) Co-Ol(A) Co-02(A) 



295 


1.874(7) 


1.925(8) 


1.993(8) 


250 


1.877(6) 


1.9357(9) 


1.9792(6) 


200 


1.892(7) 


1.9326(9) 


1.962(7) 


160 


1.889(6) 


1.9329(9) 


1.962(6) 


140 


1.882(6) 


1.9332(9) 


1.972(6) 


120 


1.879(6) 


1.9325(9) 


1.971(5) 


90 


1.918(6) 


1.9235(8) 


1.934(6) 
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FIG. 2: Left: Temperature dependence of the splitting of (4 
0) reflection|2fi| in 9 — 26 geometry (below 65K the line width 
is resolution limited). Right: The JT distortion parameter 
versus temperature as determined from the single crystal re- 
finements. The lines are guides for the eye. Inset: Profile of 
(4 0) reflection at 45K and 120K. The profile at 120K can 
be accurately fitted by two Lorentzians with the line width of 
the 45K data. 



ysis is of crucial importance for detecting the lowering 
of the symmetry. However, analyzing the single crystal 
data turns to be difficult due to pseudo-merohedral twin- 
ning and only yield proper crystallographic results when 
appropriate symmetry and twinning are used. No accu- 
rate lattice parameters can be derived from single crystal 
X-ray diffraction as they are averaged by the twinning. 
Nevertheless, the rotation and JT distortion can be ac- 
curately determined from a refinement. Once the correct 
space group is assigned, the lattice parameters can be 
more precisely determined from a powder X-ray diffrac- 
tion measurement. 

We derive three unequal Co-0 bond lengths over the 
studied temperature range: one short, one long and one 
medium. The long and short Co-0 distances correspond 
with the bonds in the ab plane, while the medium Co-0 
distances are the out of plane bonds. These results indi- 
cate a clear coherent Jahn- Teller effect associated with a 
Q2 type of distortion The right side of fig. 2 shows 
the temperature dependence of the Jahn- Teller parame- 
ter, defined as the difference between the long and short 
bond lengths normalized by their average. The magni- 
tude of the JT parameter increases with temperature, 
reaching ~ 6% at room temperature. This value is some- 
what smaller than the JT distortion of ^ 10% found in 
LaMnOs at room temperature j2^. 

The distortion of the octahedra lifts the degeneracy of 
the eg orbitals and XAS data show that the - 
like orbitals are thermally populated. The - like or- 
bitals have an antiferrodistortive ordering along all three 
directions suggesting a ferromagnetic coupling between 



magnetic sites . Such a ferromagnetic exchange inter- 
action is also supported by polarized neutron measure- 
ments 0. An antiferrodistortive orbital ordering giv- 
ing rise to ferromagnetic interactions was also found in 
YTiOs, where the t2g orbitals of Ti^~^ cations are occu- 
pied by a single electron js^]. However, for this system 
the JT distortion involves t2g orbitals for which the mag- 
nitude of the distortion is roughly five times smaller than 
for Bg-based JT systems. 

The discovery of a coherent JT distortion in LaCoOz 
offers novel insight and is unexpected considering the 
scrutiny with which this compound was studied, ft is 
clear that the JT distortion may have gone unnoticed 
in single crystal X-ray diffraction experiments because of 
twinning. Moreover, the JT distortion is not easily de- 
tected by powder X-ray diffraction because of the very 
small monoclinic distortion. 

Our measurements are in agreement with the recent 
magnetic and XAS measurements [20I Eoj . which were in- 
terpreted with a thermal activation of the IS spin state. 
Considering the fact that LaCoO^ is close to a metallic 
ground state, the transition may involve more complex 
behavior. It is observed in the manganites0, |^ that 
phases coexist of metallic and insulating states arises 
from a competition between JT-distorted and non-JT 
distorted regions in the sample. The propagation of 
strain in crystal may result in an appearance that closely 
resembles a transition. Such phase coexistence requires 
detailed real space measurements. 

We conclude that LaCoO^ exhibits a monoclinic dis- 
torted structure in the temperature interval 20K-300K. 
The monoclinic distortion state is brought about by a co- 
operative Jahn- Teller effect which triggers the long range 
orbital ordering of the Cg orbitals. The gradual increasing 
with temperature of the Jahn- Teller distortion suggests 
the thermal population of the IS state of the Co^^ ions. 
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